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SUMMARY 


This report was prepared under NASA contract NASI-11923. The scope of the 
work includes a critical evaluation of existing optical remote 
sensors for the specific problem of HCl vapour detection in the 
plumes of solid propellant rockets. 

The work performed includes, a selection of the P branch of the fundamental 
vibration-rotation band as the most promising spectral feature to sense, 
a computation of transmittance for HCl vapour, an estimation of interferent 
spectra, the application of these spectra to computer modelled remote 
sensors, and a trade-off study for instrument recommendation . 



Section 1 


INTRODUCTION 

Solid propellant rockets are expected to inject a variety of 

gases into the atmosphere. These vapours include such potentially dangerous 
gases as CO, Nitrogen oxides and HC1. The total amounts of these materials 
injected into the atmosphere by launches of different vehicles has been 
computed and is shown in Figure 1. This figure is taken from NASA draft 
report EIS-AA-724988D . ( 1972 ) . The quantities in NASA report EIS-AA- 

72-4939F agree reasonably well with those given in Figure 1. 

These gases are ejected from the rocket engine at velocities in the region 
of 8000 feet/second and at a stagnation temperature of 5,500°F. (NASA 
report N72-22790) . This hot cloud rises rather quickly due to its bouyancy, 
and through mixing and adiabatic' expansion will reach an altitude and 
temperature where it is in equilibrium with its surroundings. This 
relaxation process should occur rather quickly. Some experimental measurements by 
Forbes et al (NASA report CR-2049) /who flew an aircraft through the exhaust 
clouds of rockets , observed that the cloud temperature approached that of 
its surroundings within a few minutes of launch. Further, Block Engineering 
(USAF , Giffiths AFB report SIVS/R-23) observed that the cloud temperature dropped 
so quickly that it was difficult in their experiment to record emission from 
the cloud after the first minute. 

Following this relaxation, where the cloud comes to an equilibrium bouyancy 
with its surroundings, the usual diffusive and convective transport mechanisms 
will begin to disperse the cloud. The dispersion pattern from such clouds 
has been computed (NASA report EIS-AA-72-4988D) and a reproduction of the 
estimated HC1 concentration as a function of downwind distance is given in 
Figure 2 for the case of a low level Titan 111C abort. 

The ultimate fate of HC1 vapour is not known. However, because HC1 is a 
reactive gas, extremely hygroscopic, and highly soluble in water it is very 
probable that HCl will be removed by the atmospheric aerosol and by 
soil or water at the earth's surface. Remote sensing instrumentation could 
be used to evaluate such possible sinks, in addition to determining the 
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QUANTITIES OF POTENTIAL POLLUTANTS 
EMITTED INTO SELECTED ATMOSPHERIC LAYERS 


Note: To convert to pound*, multiply kilograms by 2.20 
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DOWNWIND DISTANCE FROM LAUNCH PAD, 
METERS 


FIGURE 2. ESTIMATED PEAK HCf CONCENTRATIONS 
DOWNWIND OF ON- PAD ABORTS. BANDS 
FOR EACH VEHICLE INCLUDE RESULTS 
FROM ANALYSIS FOR THREE ATMOSPHERIC 
STABILITY CRITERIA AND THREE INVERSION 
HEIGHT. 
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degree of pollution due to the launch of these vehicles. 

HC1 may be detected by a variety of methods including gas chromatography, chemical 
techniques and perhaps tuneable laser systems. Remote sensors which are 
capable of detecting the total burden of vapour in the sensor field of view 
without the need of active sources are however very attractive for mapping 
plumes, where large area coverage in a short time is required. 

Barringer Research Limited has developed a variety of optical correlation 
sensors which are capable of rapid and specific detection of chosen gases. The 
first such instrument developed was the correlation spectrometer (COSPEC) 
described initially by Barringer and Schock (1966) and subsequently in several 
reports and also in a recent paper by G. S. Newcomb and M. M. Millan (1970). 

This instrument is being sold commercially primarily to users interested 
in total SO^/NO^ burdens. The correlation interferometer has been under 
development for several years , and has also recently been described by Dick 
and Levy (1970) , and by Grenda et al (1971) . 

The final Barringer correlation device is a gas cell non-dispersive analyzer, 
with recently devised compensation techniques which cause the instrument 
to be insensitive to operating environmental conditions. This instrument has 
not been described in the open literature and a breadboard version is being 
fabricated at Barringer Research at the time of writing. 

In addition to this instrument there are several additional manufacturers of 
interferometers including Block Engineering and others. These are reviewed 
by Bell (1972). Furthermore, General Dynamics have over the past several 
years developed a non-dispersive infrared gas cell analyzer called a Gas 
Filter Correlation sensor which has been described by Ludwig et al (1968, 

1969) and Bartles et al (1971).. 

The scope of the present work is to evaluate the feasibility of applying 
these optical correlation remote sensing devices to the problem of 
detection of HC1 vapour in the plume of solid propellant rocket exhausts. 
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Section 2 


SPECTRA AND ATMOSPHERIC EFFECTS 


2.0 THE HC1 SPECTRUM 

The HC1 spectrum contains electronic transitions which occur in the 130o£ 
region. These are too far into the ultraviolet region of the spectrum 
to be useful for remote sensing. This is because of strong absorption 
due to 0 2 and other atmospheric constituents which remove the solar 
energy in the higher atmosphere. 

The microwave spectral region contains the pure rotational spectrum of 
HC1 , however at the pressure of the atmosphere near the earth's surface 
it is unlikely that specific detection is possible here. 

There remains the vibrational-rotational spectral structure which for the 
fundamental transition occurs at about 3.5 microns. The HC1 vibrational- 
rotational structure is shown in Figure 3. It is this spectral region 
which has been selected as the most promising for the remote detection 
of HC1 in the atmosphere. 

The HC1 transmittance spectra were computed by a local university using 
a radiative transport computer programme (Cann and Nicholls, 1972). 

The HC1 transmittance spectrum was computed for a 295°K isothermal 
atmosphere, for a 2 kilometer path of HC1 at concentrations of 0.5 parts 
per billion (ppb) , 5 ppb and 50 ppb so that the integrated HC1 
burdens were 1 ppm-m, 10 ppm-m, and 100 ppm-m. These amounts are well 
below the anticipated initial burden of more than 1000 ppm-m based 
on measurements made by Block Engineering (USAF, Griffiths AFB report SIVS/ 
R-23) . 

A sample of the HC1 computed transmittance (1 ppm-m) is shown in Figure 4. 
The transmittance due to these burdens of HC1 were used to evaluate the 
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change in spectral radiance at the sensor, and this coupled with the transfer 
function of the instrument enabled the evaluation of instrument sensitivity. 

2.1 SPECTRAL INTERFERENTS 

Two sources of spectral interferents are the natural species normally found 
in the atmosphere, and those species introduced by the solid propellant rockets 
at the same time as the HC1 target vapour. An indication of the naturally 
occuring spectral interferents was estimated from Migeottes Atlas. From this 
atlas it is evident that the R branch of the HC1 band is rather badly 
interferred with by methane, and the P branch, although somewhat "cleaner", 
has some N 2 0 interference and some water absorption lines are also present. 

In addition the atmospheric aerosol and the Rayleigh wings of distant bands 
tend to cause a general decrease in transmittance with increasing zenith distance. 
The aerosol extinction coefficients given by McClatchey et al (1972) show 
that horizontal geometry will result in a rather low transmittance. 

The result of this is a rather "warm" black body radiator against which 
the emission of HC1 will be difficult to detect. 

The interferents produced by the launch vehicle will most probably include 
aerosol material composed primarily of Al 2 C> 3 t * ie ^ to m * cron diameter 

region. In addition, surface dust will be raised by the high velocity effluent, 
and finally, it is possible that minor amounts of N 2 0 are produced. It is 
clearly a requirement of a successful HC1 detector that it be able to 
discriminate against these interferents. 

Another possible class of interferent vapours are the hydrocarbon CH strech band 
features near 3.5 microns. The degree to which these present possible 
interferent responses are difficult to predict unless each species can be 
individually assessed. Again experimentation would appear the most effective 
approach . 

So as to incorporate spectral interferents in a somewhat more realistic 
fashion, the atmospheric transmittances TR^ were obtained from Migeottes 
atlas, and this was used to estimate the interferent spectrum. This was 
chosen over the rather unattractive alternative of incorporating all of the 
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possible interferent gases (CH^, N^O, H^O, etc ^ and aerosol effects. To 
evaluate the effect of changing amounts of interferents in the atmosphere, 
the transmittance spectrum of Migeotte was modified. The simple procedure 
of taking the square and the cube roots of TR^ were chosen as the interferent 
changes. These changes were presented to each sensor model to evaluate the 
change in response caused by this interference. 

An example of a section of Migeottes atlas transmittance is given in 
Figure 5. 

2.2 ATMOSPHERIC PHYSICS OF REMOTE SENSING 


2.2.1 Measurements Made in the Downward Looking Mode 


At wavelengths below about 3y the solar radiance reflected from the earth's 
surface is greater than the thermal emission from the surface. The 
wavelength at which the reflected radiance equals the emitted radiance 
depends on the sun zenith angle, cloud cover, the surface reflectivity and 
emissivity, and on the surface temperature. In general however, this cross- 
over region is about 3 microns. The P branch centered at about 3.6y 
would normally be expected to have about 10 times more radiance emitted from 
the earth than the reflected solar radiance (Wolfe, 1965). For early 
morning, late evening, or night time launches, the solar contribution can 
be neglected, and this should be the case even during daytime launches. Certainly 
for detection requiring no more than 10% accuracy it should be possible 
to neglect the solar radiance term. 

Appendix A suggests that under certain assumed conditions the spectral 
radiance N V (S) at the instrument sensor can be expressed approximately as 

N V (S) = e E B v (T e ) Tr (S) + (1 - Tr y (S ) ) B v (T ft ) (1) 

where e is emissivity, B the black body radiance at temperature T and T (S) 

v rv 

is the spectral transmittance through the atmosphere to S. 




| WAVENUMBER (cm -1 ) 


ATMOSPHERIC TRANSMITTANCE FROM MIGEOTTES ATLAS 


FIGURE 5 
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The subscripts E and A refer to the earth's surface and the atmosphere 

respectivelv . it is evident that if e B (T ) = B (T ) one may write, 

E u E v A 

N (S ) . = B (T ). 

V V A 

which shows that the spectral radiance is independent of the atmospheric 
transmittance and gas detection is impossible with, these assumed conditions. 

Figure 6 shows that typically each morning and evening there is a time when the 
earth's surface temperature equals the air temperature. If the surface 
emissivity is approximately 1, there will be times during which normal 
downward detection of atmospheric vapour is not possible. However, 
measurements even at these times may be possible if one can harness naturally 
occurring variations in the albedo. This so called "ground chopping" 
technique has been described (Zwick and Millan, 1971) . 

Ground chopping is a technique which relies on the spatial variations in 
surface radiance (caused by emissivity and temperature variations) . If the 
sensor "footprint" on the surface is small there will be a time varying 
radiance from the surface as the footprint spatially scans the non-homogeneous 
surface. If the sensor electronics is band limited so that only the time 
varying spectral radiance, n , is accepted (AC coupled) then the spectral 
radiance at the sensor becomes 

n (S) = n (o) Tr (S) (2) 

V V V 

The emission from the atmosphere will only vary slowly in time as the spatial 
extent of the pollution cloud is large with respect to the footprint size. 

Such slow variations may be rejected by the electronics. Only rapid 
atmospheric temperature variations could result in a radiance fluctuation which 
could be passed through the filter, but these are expected to be small. 

Ground chopping has not been demonstrated in the infrared but it seems worth 
consideration and testing for many terrain types. 
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Downward looking sensors can however be used for HC1 detection when the 
atmosphere and the earth's surface differ in radiance, as is normally expected 
near mid-day or near mid-night. 

2.2.2 Measurements in the Upward Looking Mode 

An upward looking system has the advantage that the rather small radiance 
of space can be used as a background against which the emission of HC1 
at atmospheric temperatures can be sensed. The other upward looking 
geometry is to use the sun as a source, however for plume mapping the 
additional complexity of maintaining solar sighting is undesirable. 

For the case of an upward looking sensor not looking at the sun, the 
spectral radiance can then be written as 

N (S) = (1 - TR y (S)) B^ (T ) (3) 

If the transmittance is expressed in terms of the target vapour 
transmittance times the interference vapour plus the aerosol transmittance, 
one obtaines 

N y (s) = (1 - TRj(S) TR^(S)) (T ) (4) 

i 

In the upward looking mode the interference transmittance is that integrated 
along the line of sight of the sensor through the whole atmosphere. 

Obviously for angles approaching the horizontal, TR*. becomes small. 

A rather crude extrapolation of Etermans (1970) data to 3.6 microns to 
the case of a visual range of 10 km suggests an extinction optical thickness 
of about 0.2. If one makes the crude assumptions that the atmosphere 
is isothermal up to 1 km, and so cold above that the radiance from greater 
altitudes may be neglected, then the effective extinction thickness is that 
between 0 and 1 km, or about 0.1. Further, according to Hodges (1972) , 
about 10% of the extinction at 3.6y is due to absorption (and thus 
presumably can be re-emitted) , so that the effective optical thickness due 
to this 1 km isothermal absorber would be about 0,01. 
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Although this thickness should increase with angle, according to sec 0 
approximately, rather large zenith angles should be possible before 
detection is severely limited. These calculations are very crude, 
and allowance is not made for the aerosol generated by the launch vehicle. 
These effects will require experimental evaluation, but the indications 
are that this operational mode should also be feasible -provided the sensor 
has sufficient sensitivity. 
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Section 3 


OPTICAL CORRELATION REMOTE SENSORS 

3.1 INTRODUCTION 

The theory and operation of some remote sensors are outlined in this 
section. Since a comparative study is to be made, several features 
have been chosen as indicative of sensor capability. These features include 
sensitivity, specificity, observation and integration time, size— weigh t-and- 
complexity, and cost. As these are features common to each instrument, 
some discussion of these will preceed the actual instrument discussion. 

3.1.1 Sensor Sensitivity 


The sensor whose sensitivity to HC1 is greatest is the most desirable. 
Sensitivity to an integrated burden of approximately 30 ppm x 100 meters 
is about the upper limit expected in the ground cloud, and to be useful the 
sensor should be capable of sensing HC1 burden several orders of magnitude 
lower than this. 

The instrument sensitivity can be specified in terms of the power change 
at the detector caused by the presence of the HCl vapour. 


The power P, reaching the detector when no HCl vapour is in the atmosphere is 


P 


= / N v (S) xAft dv 


(5) 


where N^ (s) is the spectral radiance at the sensor collector when no HCl is 
in the atmosphere 

A is the collector area 

ft is the solid angle of acceptance of the collector - 

x is the spectral efficiency of the sensor through to the detector and 

for convenience is written as the product of the average transmittance 

x times the relative spectral transmittance x . 
o v 
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The function x is the relative instrumental transfer function, 
v 


The change in power, S, at the detector caused by the presence of the 
HC1 vapour is 


S 


= iKa s) ~ N v CS)) 


T 

V 


Aftx dv 
o 


v 

Where N^CS) is the spectral radiance at the sensor collector when HC1 
vapour is present. 


(6) 


If one assumes a detector noise limited device, as may be expected for an infrared 
sensor system, the noise power level, NEP is given by 


NEP 


( A d Af) 


1/2 


D* 


(7) 


Where A, is the detector area 
d 


Af is the effective bandwidth which is related to the integration time Tv 

m 

by T m = l/(27rAf) 


and D* is the figure of merit of the detector. Hence the signal-to— noise 
may be written as 


S/N 



x dv 

V 


] 


Afix D* / (A.Af 1 
o ct 


1/2 


( 8 ) 


For the sensor types to be considered, the variables are A,.fi, x_, and x and A. . 

° I v d 

The same detector type will be assumed for each instrument and the integration 
time will be cassumedi. equal in each case. 


3.1.2 Sensor Specificity 

The specificity is defined as the degree to which the sensor response is 
insensitive to interferent species. In the case of a classical single entrance 
slit, single exit slit spectrometer the specificity could be rather simply related 
to the instrument resolution. The higher the resolution the better the 
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rejection of interferent species, up until instrumental bandwidth is narrower 
than .the target spectral feature. This same concept is true of correlation 
instruments. However, with greater specificity goes a lesser sensitivity. 

In the case just cited, the narrower spectrometer slits have less radiant 

throughput or etendue. In fact the resolution-etendue product 

is a constant for grating spectrometers (.James and Steinberg, 1969) . 

In the case of optical correlation instrumentation the effect of reducing 
resolution so as to increase sensitivity results in a loss of interferent 
rejection which is not so obvious to evaluate as with conventional 
devices. It is possible by judicious choosing to accentuate certain tarcjet 

I 

spectral features which have few interferents , and to reduce the influence 
of those regions which have serious interferents. 

It is because of this complex spectral trade-off that computer programmes 
were devised to simulate the response of each sensor and to evaluate the : 
ability of each to reject a simulated spectral change in atmospheric j 
radiance. i 

! 

3.1.3 observation and Integration Time 1 

! 

Commercial instrumentation is usually designed so as to enable reconstruction 
of the spectral radiance of the scene being viewed. Hence it is usual for 
spectrometers and interferometers to require of the order of seconds to jcomplete one 
spectral scan (interferogram in the case of interferometers) . Bell (197(2) gives 
a review of commercial interferometers available in about 1970. 

With correlation instrumentation however the situation is quite different. 

Because the instruments are by nature multiplexed (observe all desired j 
spectral features simultaneously) it is only necessary to oscillate 
between two pre-chosen multiplex positions or regions in order to ! 

execute a measurement. This oscillation can be performed in abo ut one 
second with correlation interferometer sensors, about 10 milliseconds w^th 
some correlation spectrometers, and several thousandths of a second for Isome 
correlation gas cell analyzers. 
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It is this possibility of very fast observational measurement which makes 
the optical correlation instruments ideal for airborne sensing, or 
for studying the dynamics of flow systems such, as cloud dispersion where 
fast large area coverage is required. 

3.1.4 Size, Weight and Complexity 

Other trade-off parameters to consider in evaluating an instrument for a specific 
task include size, weight, complexitv and ruggedness. The sensor best suited 
for plume dynamic coverage must be rugged and portable. Again optical correlation 
sensors score heavily in these categories. 

3.1.5 Cost 

Obviously the most cost effective instrument is that one which performs 

the desired measurement at the least cost. A cost evaluation for each sensor 

will be given in a following section. 

3 . 2 CORRELATION SPECTROMETER 

3.2.1 Introduction 


Several optical correlation sensors have been developed which offer the 
possibility of sensitive, rapid remote monitoring of integrated gas burdens 
along the line of sight of the sensor. Correlation spectroscopy instrumentation, 
using grating spectrometers .have been discussed by Barringer and Schock (1966) , 
Ludwig et al C1968) , Williams and Kolitz (1968), Millan et al (1969) j Davies 
(1970) , Newcomb and Millan (1970) , Millan (1971) , and Moffat et al (1971) . 

The closely related technique of Hadamard Spectroscopy has been discussed by 
Ibbet et al (19681, Decker and Rarivet (1968) , Nelson and Fredman (1970). ! The 

I 

use of multiple entrance slits for added throughput or etendue, has resulted 
in a coming together of classical spectroscopy and Hadamard Spectroscopy. ' 
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The name of correlation spectroscopy arises from the fact that the observable 
chosen as output is the variation of the correlation function between a 
mask, or array of exit slits, situated in the exit plane of a spectrometer, 
and the dispersed power spectrum, produced in the spectrometer. The correlation 
function between the dispersed power spectrum of the incoming light and the array 
of slits is the power in watts passing through the mask, from the dispersed 
power spectrum to a photo-sensitive device situated behind the mask. 

The correlating function may be of the matched filter type to match the absorption 
spectrum of the absorbing target gas, or the filter may be optimized to give 
a minimum rms response due to anticipated inter ferents. A correlation 
spectrometer is designed to detect the quantity of a specific gas present in 
the line of sight by displaying the dispersed spectrum of incoming radiation 
at the exit mask. The mask structure allows alignment to specific absorption 
bands of the gas under investigation. 

Modulation is produced either by oscillating the spectral position relative to the 
correlation mask, or alternatively by cyclically moving the mask position as seen 
by the detector. The modulation is detected by a detector and electronically 
processed to provide an output proportional to the concentration pathlength 
product. The optical layout of one such device is shown schematically in Figure 7 

3.2.2 Operational Description 

The dispersive correlation spectrometer consists of fore-optics, a modulator, 
grating spectrum disperser, exit mask, detector and associated electronics. 

The modulator may consist of a pair of refractor plates mounted in the end of one 
of the arms of a tuning fork, which is driven at its resonant frequency (normally 
in the 100 - 1000 Hz ranged. The movement of the plates causes the light beam 
to "jump" cyclically between two predetermined positions, (1) and (2). The 
distance between the positions: is determined by the refractor plates' 
thickness, angle between plates and wavelength of radiation. (In another 
configuration the masks are mounted on a rotatngg disc which alternately passes a 
pre-chosen set of spectral wavelengths to the detector) . 
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SCHEMATIC LAYOUT OF A CORRELATION SPECTROMETEF 
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The dispersed spectrum in the exit plane of the polychroma tor is a multi- 
wavelength image of the entrance slit. The position of any predetermined wave- 
length with respect to a set of fixed axis in the exit plane is governed 
by the grating angular position. 


The mask or array of exit slits is situated in the exit plane in front of 
a detector. The response of the detector-preamps, is a variable quantity. In 
order to control the amount of current from the detector under varying light 
conditions, an automatic gain control (AGC) is normally- used, which together with the 
synchronous detector perform the following functions; 

(a) One of the positions of the spectrum with respect to the exit mask and with 
incident power P is selected to control the AGC. 

(b) The AGC loop is phase-locked to this position and a fixed signal is 
selected according to the expected use of the instrument. The gain 

is automatically adjusted so that the response to mask position 1 is always 



(c) The second position of the spectrum with respect to the mask with incident 
power P 2 generates a detector signal S 2 - 

(d) The process is repeated for each oscillation of the fork. 

(e) A reference signal is obtained directly from the fork oscillation and a 
gate waveform is derived to be used in the synchronous detector. 


(f) S 2 will vary with, respect to as a function of the characteristics of the 
dispersed spectrum. 

(,g) The product detector waveform is passed through an integrator and a signal 
output is presented to the recording device. 


-•14 - 


21 



The difference signal is given by 
S -(Sj - S 2 )d 

V 

where (fixed) a G a — 

P 2 

and S 2 aP 2 G aS i “ 


The signal current then becomes: 



where d is the duty cycle of the gate waveform and will be included into 
x the optical efficiency in the remainder of this report. 

The AGC circuit is provided for two main reasons: 

(a) to provide automatic normalization between output current in mask positions 
(1) and (2) , thus compensating for multiplicative changes in radiance level 

(b) to maintain constant the output current of the detector from one of the 
(mask-spectrum) positions thus providing a normalizing reference and 
easing the power supply demand on the detector. 

The AGC then provides a feedback loop which adjusts the effective gain. 

The Ctendue achievable with a dispersive spectrometer is limited by the 
entrance slit width in the direction of dispersion, due to whatever resolution 
is required. Slit heights can, on the other hand, be made fairly large. 
Efficiency in transferring gas modulated radiation to the detector is 
reasonably good. Many of the optical surfaces may be anti-reflection coated; 
the main loss is in the grating efficiency. The use of a mask in the 
plane of the spectrum allows the best possible modulation to be obtained, 
within the limits of resolution and the requirement for avoiding regions 
subject to interferences. 
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Attenuation of interferent takes the form of avoiding affected spectral 
regions where possible, and balancing them against each other where not. 


Hence the device transfer function is 


V F V M V1 and t v 2 = f v m v 2 . 


The device response and S 2 when only interferents are present are 

. S 1 = g.:-.a.q.t 0 (S) . dvj 


S — G . A.fi.XO 


[ Jk (S) * T v2 dV ] 


(10 


F v is the interference filter function which isolates the P branch of HC1, 
and M v 2 are the two mask functions obtained from the convolution of the 
entrance slit and the exit slits. 


The response S - S* - S* will be an offset response; constant to the extent 
the instrument is insensitive to interferents. 


The downward looking HC1 response (from equation 1) can be written 


as S. = s g - 


S d = GAOTq 


J"v p v [wV Tr i (S) TR y CS) 


+ B v (T a ) (1 - TrJ(S) Tr^ (S))J dv 


jVv ( e E \ ( V Tr J < S > Tr y (S) + B v (T a ) U 

V2 


- S' 


Tr^(S) x Tr^(S) dv 
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This reduces to 



GA^t o 


r 

M.F 
v V 

V 1 


d 

X dv 
v 


f M F X d dv 
J V V V 

v 2 


( 11 ) 


where 


X d = e B (T ) Tr 1 (S) (1 - Tr g (S)) - B (T )Tr X (S) (1 - Tr g (S)) 

vEvEv V vAv V 


When ground chopping is employed the above equation has a new value of X^ 
given by 


x gc 

v 


i 

c Tr 

v 


(S) (1 - Tr g (S) ) 
v 


( 12 ) 


where C is the a.c. component of £ E B V • For th e computation C was chosen 

to be 10% of e B (T ) . 

E V E 

In the upward looking mode the signal response due to HC1 vapour is obtained 
from equation 11, but where X^ is (using equation 4) 

X U = B (T ) Tr^S) (1 - Tr g (S)) ^ 13) 

V V A V V 

3.2.3 Correlation Spectrometer Parameters for HC1 Detection 

The masks were chosen so as to transmit six H d Cl 37 P branch, lines (P 2 through 
P 7 ) , and the spectral width of the entrance and exit mask were chosen 
to be approximately 2 cm' 1 wide (25$) . Mask 1 was placed so as to centre the 
transmission function on the line location, and mask 2 was chosen to be centered 
3.5 cm 1 toward larger wavenumbers. 
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Other instrumentation parameters were 


-3 2 

Aft = 1.4 x 10 cm Sr 

T = 0.13 

° 2 
A , = .16 cm 

d 1/2 -1 
D* = 10" cm Hz W 

Af = 1 Hz 

The- computer programs are discussed briefly in Appendix D. The results of 
this modelling are given in the next section. 

3.3 CORRELATION INTERFEROMETER 

3.3.1 Introduction 

An interferometer transfers to the detector not the radiation at 

a given wavelength as does a dispersive spectrometer, but rather the interferogram 

Cfourier transform) of the radiation. The relationship between the 

two is outlined in Appendix B. The output takes the form of cyclic variation 

of the amount of radiation transmitted, going through approximately one 

cycle or "fringe" for each wavelength change in interferometer delay. 

Ordinarly interferometry has advantages over dispersive spectrometry in 
two respects . 

The radiance throughput or etendue, (.the amount of available radiation which 
the device can accept) is usually greater for an interferometer than for a 
dispersive spectrometer. The interferometer is also completely multiplexed 
compared to the correlation spectrometer which often uses only the main spectral 
features . A Hadamard spectrometer in principle has the multiplex and etendue 
advantages of ordinary interferometers. Block Engineering have on the market 
a family of ordinary Michelson interferometers. 


18 


25 



There are however certain forms of interferometers which, enjoy additional 
advantages over ordinary michelson or Fabry-Perot interferometer designs. 
Connes (J.958) has devised a spherical Fabry-Perot design which at high 
resolution has an etendue advantage. Hansen (JL941) and subsequently many 
others have shown that a Michelson interferometer may be field widened so 
as to have an additional radiance throughput advantage over ordinary 
Michelson interferomters . Reviews of this principle may be found in Steel 
C1967) , James and Steinberg (1969) . . 

Simply stated, field widening occurs when there is an uncompensated thickness 
of refractive material in one arm of the interferometer. The field 
widened condition occurs when the image position of one mirror is coincident 
with the second mirror position but the additional refractive material 
in one of the arms results in a non-zero delay. The interferometer 
depth-of -focus for a plane parallel slab of refractive material in one arm 
has been computed by Hilliard and Shepherd (1964) and derived by Zwick 
and Shepherd (1971) . Barringer Research Limited has developed a correlation 
interferometer (Grenda et al, 1971) which scans a small delay region centred 
about the field-widened position. In this case the scanning is achieved 
by rotating the refractor plate. In the case of the Barringer instrument the 
usual spherical aberration and defocussing of a wide angle interferometer 
are present, and in addition there is astigmatism. However the use of high 
refractive index slabs results in acceptable fringe visibilities over a few 
millimeters of delay scanning distance (Zwick et al, 1971). As this scan 
region can be centred at any pre-chosen delay position it is possible to scan 
through sufficient independent spectral elements to perform a spectral 
analysis on the interferogram segment and so to achieve good specificity 
to supplement the high sensitivity of the field-widened device. 

Correlation interferometry is a technique for extracting specific desired 
information from the radiation interferogram as efficiently as possible, 
in order to keep to a minimum uncertainties due to all sources. In this case 
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correlation may be done, on interferograms. obtained with, conventional Fourier 
instruments such, as the Block, spectrometer , or a set of pre-determined 
weights (correlator! may be used which are applied to the interferograms in real 
time with rather simple hardware. 

3.3.2 Optical Arrangement 

The basic configuration for a correlation interferometer is shown in Figure 8. 
Operation of the instrument consists of varying the interferometer delay by 
either moving one of the end mirrors or by rotating the scanning refractor 
plate. The latter technique offers great simplicity in comparison to 
scanning of one of the end mirrors, as there is no possibility of spoiling 
the alignment by tilt effects. The scan is monitored by a reference 
interferogram and by means of an encoder fixed to the scanning plate. Ordinary 
interferometers scan through zero delay whereas the refractively scanned 
device merely uses the encoder pulse to sense the beginning of the scan. 

A filter is used to isolate the spectral region of interest. If several 
separate regions are of interest, multiple detectors may be used to advantage 
in conjunction with dichroic beamsplitters and spectral filters. 

3.3.3 Signal Processing 

Processing of the resultant signal is carried out in two stages. In the first 
stage, fringes from the reference source are used to heterodyne the high frequency 
carrier present in the signal while at the same time automatically compensating 
for small irregularities in the scan or drifts in the optical paths. 

In correlation interferometry a very high degree of specificity is obtained from the 
resulting preprocessed signal by the application of a linear digital filter. This 
filter consists of a set of weights chosen in such a manner as to provide an output 
which measures quantitatively the presence of a specific spectral signature, 
but which is insensitive to the interfering spectral signatures. This filter 
may simultaneously be optimized for the rejection of various types of synchronous 
and non-synchronous noises . Appendix B outlines the relationship between the 
the spectra and the interferograms, while Appendix C describes the mathematical 
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Figure 8. OPTICAL LAYOUT OF AN INTERFEROMETER 
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techniques for selecting the weighting functions. In setting up 
for operation, sample interf erograms are observed for as- wide 

a range of interf erent conditions as is considered probable under operation 
conditions. The techniques of Appendix C are then used to determine the 
weighting function (s). 

One powerful feature of the technique is its capability for economically 
measuring several different signatures whenever these provide information 
over overlapping regions of the inter ferogram, simply by the simultaneous 
application of different weighting functions. In cases where the signatures 
occur in different spectral regions, separation prior to the detector as 
mentioned above eliminates crosstalk between outputs; however where 
crosstalk is not a problem the radiation may be received by a single 
detector. 

A somewhat emperical approach will be used to express the S/N response of an 
interferometer sensor for HCl detection. The methods used to evaluate noise in 
Fourier spectroscopy are reviewed in James and Steinberg (1969). Basically the 
difficulty is that the detector receives all of the radiant energy all of the time 
(true multiplex 'advantage) but the time varying component is only that of the 
Fourier transform of the spectral radiance at the given delay and scan rate. The 
detector noise is fixed and easily evaluated but noise due to radiance 
changes (multiplex disadvantage) are indeed difficult to express. 

A signal response is therefore here evaluated as an average value over the 
interf erogram correlation region. A modulation efficiency, m, is derived 
from computed interferograms , and is defined as 



Where I (Xl is the interf erogram value at the delay region X, < > 
denotes average value, and g and i refer to the interf erograms when HC1 
vapour is included and when it is not included respectively. Values of 
m were evaluated directly from the. interf erograms. 


Once these modulation efficiences are defined one can write the response 
due to HC1 vapour by using equations 6 and 14, and including the interference 
filter , 


S 


AfitQm 


(N^(S) - N g (S) ) . F dv 


10 



(15) 


Values of N (S) and N g (S) for the downward looking mode are obtained readily 
v v 

from equation 1 and in the case of interferents only TR^(S) = TR^(S) 

where as for the case including HC1 the transmittance becomes TR^(S) = TR^(S). 

TR g (S) . 
v 


In the downward mode, but using ground chopping the response is 

r 


qc 

S = m.A(2x f 


|C Tt (S) (1 - TR 9 (S) ) . F dv 

R V v v 


(16) 


Where, as before, C was set equal to 10% of e B (T ) , the earth's average 

.fcj v 

spectral radiance. 

Finally the upward looking response is (from equations 4 and 15) 


S = mAfix 0 


B (T ) TR 1 (S) (1 - TR 9 (S)) F .dv 

v A v v v 


(17) 
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These response signals- were computed for various HC1 burdens, and 
for differing interferent transmittances. In addition the response was 
compared with, response due to detector noise as. given by equation 7. 

The results of the calculated response signals is given in the next section, 
and the programmes used are described in Appendix D. Some work was done to 
evaluate the best region of the interf erogram to perform the HC1 detection. 

The criteria chosen was to find the region which gave the optimum 
sensitivity but at the same time, the best specificity. It was found 
that scanning the .065 to 0.1 mm delay region gave the best results. 

The choice of region was not especially critical so that slight changes 
in the choice of operating region did not severely change the results 
(provided that weighting functions were derived for that region) . Approximately 
equal sensitivity would be achieved without the weighting function 
correlation, but not nearly such good interferent rejection could be possible. 

The interferometer parameters were as follows: 

-1 2 

Afi = 1.1 x 10 cm Sr 

x = 0.06 

° 2 
A^ = 0.16 cm 

D* = 10^ cm Hz^/2 w 1 
Af = 1 Hz 

3 . 4 CORRELATION GAS CELL ANALYZER 

The non-dispersive infrared gas cell analyzer has been in use for many years. 
Luft (.19431 first constructed such a device and a succession of workers 

have since used similar systems (Martin, 1953; Smith and Pidgeon, 1963; Goody, 
1968; Ludwig et al , 1969). The basic theory of operation of these instruments 
is identical, the differences being in the optical arrangement of the 
components, the methods used for modulating the incoming signal and the 
electronic signal processing. Commercial instruments have been developed by 
General Dynamics and Barringer Research Limited, and perhaps others. 
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Figure 9 shows the bas.i.c layout of one such correlation gas cell analyzer 
system. The incident radiation is divided at the beamsplitter and 
passes through one or the other of two gas cells, Cl or C2 and onto one of 
the detectors. One of the gas cells (Cl say) contains a sample of the same type 
of gas as is to be sensed, while the other contains a spectrally inactive 
gas. The optimum amount of gas is experimentally (or computer) evaluated to 
give maximum response. (In another version of this instrument the radiation 
is alternately passed through cells Cl and C2 and onto a common detector.) 

At a wavelength where the gas absorbs strongly, radiation passing through 
cell C2 (containing no target gas) is attenuated depending on the 
quantity of target gas that the radiation has passed through 

before reaching the instrument. However when the target vapour radiation passes 
through Cl, the radiation is largely removed regardless of the amount of gas 
prior to the instrument. In this path it is much less affected by 
the quantity of target gas prior to the cell. 

A chopper is used to modulate the input source radiance to a frequency 
best suited to the detector (in the 3 k Hz region for cooled PbS) . The 
chopper also chops an additional source whose purpose is to provide a 

synchronous signal to be used in. .a product detector for synchronous 
monitoring of the source radiance at each detector. The difference in the 
output from these two synchronous source signals is an output which varies 
in a manner indicative of the amount of target vapour signature at the- 
sensor. The signal processing is shown in a block diagram form in Figure 10., 

Normally the two detector responses would vary with ambient changes in 
responsivity as the two detectors age, change temperature, etc. Further 
the response could depend to some extent on the thermal stability of the 
two arms of the instrument. To help remove these undesirable effects 
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a reference lamp is included in the system. The detector then senses both 
the source and the reference lamp. Another signal, synchronous with the 
reference lamp is used to synchronously extract the reference lamp response 
from each detector and the difference between the two detector reference lamp 
responses is held constant by a feedback, loop, which adjusts the gain of one 
detector. 

In this way the two detector responses are forced to "track" each other. 

In the case where the radiation is alternately passed through Cl and C2 to 
a common detector, the optical paths may be forced out of balance by thermal 
emission changes in the two optical paths, or dust specks etc. The two 
detector forced tracking system eliminates these effects. An additional 
feedback loop is used to keep the average reference lamp radiance equal 
to the average source radiance so that differences due to changing detector 
responses with input power are minimized. 

The correlation gas cell analyzer has a high degree of specificity since 
absorption of radiation at wavelengths other than those absorbed by the 
target gas affect both outputs equally, and so does not affect the 
difference. Rejection of interferents is very good when both target 
and interferent gases have a large ratio of line separation to line 
width as in the HC1 case. Laboratory work by Bartle et al (1971) has 
demonstrated the rejection capability of a correlation gas cell analyzer. 

The output signal from a correlation gas cell analyzer may be written 
from equation 6 : 


S = Afix 0 

V ' 

where x is the normalized transmittance through cell Cl and x the 
cl cz 


(N (S) x . - N (S) x ) dv 
v cl v c2 
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transmittance through cell C2. x^ was obtained from a radiative transport 
calculation through a cell 2 cm long containing 0.1 atmospheres of HC1 cell x 2 
was set equal to 1. The signal duty cycle (time during which source is being 
sensed) is included in the calculations. 


The values for source radiance at the receiver (si were obtained as before 
from equations 1, 2, 3 , and 4. 


The correlation gas cell analyzer parameters used in the computations 
presented in the next section were, 


2 

Aft r 0.1 cm Sr 

x = 0.1 

° 2 
A = 0.16 cm 

® 11 1/2 -1 
D* = 10 cm Hz 7 W 

Af = 1 Hz 


In evaluating the signal to noise ratio for the double element 
detector system the noise equivalent power expressed by equation 7 can 
be used, but a factor of / 2 is required to account for the rms addition 
of the two individual detector noises when the signal difference is formed. 
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Section 4 


SENSOR TRADE-OFFS FOR HC1 DETECTION 
4.1 HC1 SENSOR DETECTION SENSITIVITY 

Results of the computation outlined in sections 3.2, 3.3 and 3.4 for the three 
optical correlation instruments are given in Table I. The cases evaluated are: 

(a) upward looking geometry using a computed HC1 transmittance, and 
Migeottes spectral transmittance of the whole atmosphere. 

(b) a downward looking configuration which employs ground chopping. In this 
case only the HC1 transmittance was considered, as it is not necessary 
to look through an atmospheric amount of interferent vapour in this 

case. The ground chopping radiance is taken to be 10% of the surface average 
radiance. 

Two different integrated amounts of HC1 were evaluated; these correspond to 
a 0.5 ppb HC1 concentration over a 2 km pathlength (1 ppm-m) and a 50 ppb 
concentration over a 2 km pathlength (100 ppm-m) . 

As expected, the response is not quite linear, although at these low 
concentrations the non-linearity is rather small. 

It was beyond the scope of this work to evaluate the response at very large burdens 
of pollution. It is expected that useful measurements of integrated burdens of 
several orders of magnitude more HC1 (thousands of ppm-m) are readily achievable. 

The conclusions from the table are that the gas cell analyzer is the most 
sensitive device, followed by the interferometer and the least sensitive device is the 
correlation spectrometer. It Is possible that the geometrical parameters for these 
instruments could be changed so that the numbers would be changed somewhat. 

However, the parameters in the calculation are those of existing hardware 
devices. Changes by factors of 2 are possible, but orders of magnitude 
changes would be rather difficult. 
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Instrument 

Operating Mode 

HC1 Vapours 
and Interferents 

— — — — — — — — ——4 

Computed S/N 

(Assuming Detector Noise And 
and 1 HZ Bandwidth) 

Correlation Gas 
Cell Analyzer 

Vertically 

Upward 

1 ppm-meter HC1 
and Migeotte 
Transmittance 

_ 


— 11 - , -'l 

0.6 


Vertically 

Upward 

100 ppm-meter s HC1 
and Migeotte 
Transmittance 

53 

' 


Downward with 
Ground Chopping 

1 ppm-meters HC1 
No interferents 

0.13 


Downward with 
Ground Chopping 

' 

100 ppm-meters HC1 
No interferents 

12 

Correlation 

Interferometer 

Vertically 

Upward 

1 ppm-m HC1 and 

Migeotte 

Transmittance 

6.3 x 10~ 2 


Vertically Upward 

100 ppm-m HC1 and 
Migeotte ■ 
Transmittance 

i — 1 

« 

00 


Downward with 
Ground Chopping 

1 ppm-m HC1 no 
interferents 



Downward with 
Ground Chopping 

100 ppm-m HC1 no 
interferents 

3.6 

■ 

Correlation 

Spectrometer 

Vertically 

Upward 

1 ppm-m HC1 and 
and Migeotte 
Tr ansmi ttance 

1.5 x 10~ 2 j 


Vertically 

Upward 

100 ppm-m HC1 and 

Migeotte 

Transmittance 

1.4 


Downward with 
Ground Chopping 

1 ppm-m HC1 no 
interferents 

3.4 x 10" 3 


Downward with 
Ground Chopping 

10C ppm-m HC1 no 
interferents 

3.1 x 10 _1 


TABLE I 

HC1 SENSITIVITY OF OPTICAL CORRELATI ON INSTRUMENTATION 30 

























T he sensitivity advantage of the correlation gas cell analyzer stem mainly from 
two factors. Firstly, the average optical transmission is approximately two 
times greater than the interferometer as the latter device returns half the radiance 

back to the source. Secondly, the gas cell device modulates all wavelengths in 
phase whereas the interferometer working at an interferogram HCl peak response has some 

of the signal (especially here were isotopes exist) out of phase. 

4.2 HCl SENSOR SPECIFICITY 

The computed effect of interferent changes were performed by evaluating the 
instrument response changes which resulted when the same 1 ppm-m of HCl was 
present, but when the interferent spectrum changed. The change in interferent 
spectrum was arbitrarily chosen to be that which results in an atmospheric 
transmittance which is the square root of the Migeotte atlas transmittance. This 
should be very crudely equivalent to an air mass change of 0.5. 

This spectrum was not chosen so much because of its relationship to the 
real physical world, but rather because of ease of evaluation. Since 
this same interferent spectrum was presented to each sensor it gives some 
indication of the specificity of each sensor. The results of this 
computation are presented in Table II. 

As seen from Table II , the interferometer is by far the most specific of the 
sensors evaluated. This result is due to the fact that the weighting function 
correlator is computed so as to give the minimum offset, in the rms sense, to 
a given library set of interf erograms . If the library set of conditions brackets 
(or even better includes some of) the range of interferent changes encountered, 
then they are "anticipated" interferents and can be very effectively discriminated 
against. 

The other two sensors (spectrometer and the gas cell analyzer) have fixed correlating 
functions; a metal disc or mask in the case of the spectrometer, and an HCl gas cell 
in the other case. These fixed correlator functions are reasonably good 
at discriminating against interferent effects, but not as good as. the 
interferometer. 
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4.3 SENSOR SIZE, WEIGHT, COST AND COMPLEXITY 

A ballpark- estimated cost for fabricating each, of the three optical sensors 
for HC1 detection is. shown in Table III, along with some size, weight and 
relative complexity figures. These are estimated costs. 

4.4 THE BEST HC1 SENSOR 


As has been shown in Sections 4.1, 4.2 and 4.3 the best sensor for HC1 detection 
depends on the criterion chosen. The correlation gas cell analyzer has the 
best sensitivity, but the interferometer has the best specificity and the 
spectrometer has perhaps the best ruggedness. 

If the gas cloud can be monitored from a ground based mobile van with the 
sensor pointed in the upward direction, the interferometer would be recommended, 
great specificity of this sensor is valuable because the ground cloud plume 
may contain a great many interferent species. The time available for 
burden measurement from a low speed surface vehicle could be reasonably 
long (say 1 second) and the plume could be mapped by elevation scanning and 
by van mobility. 

The correlation gas cell analyzer has good sensitivity however, and 
seems somewhat better if the HC1 cloud is to be monitored to the lowest 
possible detectable amounts. In this case the need for high interferent 
rejection has been sacrificed for the gain in sensitivity. 

If the HC1 cloud is to be monitored by looking downward from the airborne 
platform it will either be necessary to observe when the atmosphere has 
a known temperature different from the earth's surface, or else 
possibly employ the fluctuations in source radiance from the earth’s 
surface. If the sensor footprint moves quickly so as to obtain large 
area coverage from a quickly moving platform, the fast response time 
achievable with the correlation gas cell analyzer would be useful. 

The alternative is to use image motion compensation techniques which would 
add complexity and cost to the system. 
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Instrument 

Operating Mode 

Spectral Change 

Equivalent HC1 
Response 

Correlation Gas Cell Analyzer 

Vertically Upward 

1 ppm-m HC1 with Migeotte Spectrum, 
changed to 1 ppm-m HC1 with square 
root of Migeotte Transmittance 
25 0. 5 air mass 

cs 15 ppm-m 
HC1 

Correlation Spectrometer 

Vertically Upward 

as above 

2? 20 ppm-m 
HC1 

Correlation Interferometer 

Vertically Upward 

as above* 

~ .04 ppm-m 
HC1 


The linear digital filter weighting function was in this case evaluated from Migeottes atlas transmittance 
and from cube root of Migeottes atlas transmittance. The optimum filter from this spectral set was applied 
to the case where the interferent spectrum corresponds to the square root of Migeottes atlas. 


TABLE II 

OUTPUT CHANGES DUE TO INTERFERENT EFFECTS 














1 Instrument 

j 

Size 

Weight 

Cost* 

Relative 

Complexity 

Ruggedness 

1 

Gas Cell Analyzer 

2' x 2' x 8" 
plus 3' telescope 

50 lbs 

$32,000.00 

Moderate 

Moderate 

Correlation 
. Interferometer 

18" x 12" x 12" 
plus 3' telescope 

30 lbs 

$65,000.00 

High 

Moderate 

l 

Correlation 
! Spectrometer 

2' x 1' x 8" 

30 lbs 

$30,000.00 

Moderate 

Moderate-high 


* ESTIMATE ONLY BARRINGER RESEARCH SENSORS 


i 


TABLE III 

HC1 SENSOR PARAMETERS 
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APPENDIX A 


Radiative Transport 

In the thermal infrared spectral region the transmission of radiation through the 

atmosphere is attenuated due to absorption by atmospheric gases and aerosols, 

and due to scattering by airborne particulate matter. In addition however there 

is a competing emission from the gases and aerosols which tends to offset the 

absorption radiance losses. One may express (Armstrong and Nicholls, 1972) the 

change in radiance dN (S) as follows: 

v 

dN (S) = -N CO) dx + J (S) dx 
v v v 

where dx = k^p ds is the element of optical depth, k^ is the mass 
absorption coefficient, p the density, ds the element of path length. The mass 
absorption coefficient depends on the molecular parameters, and on the 
atmospheric temperature, pressure and composition. 


J (S) is the source function and is given by J (S) =j (S}/k and j (s) is the 
v 3 v v v J v 

mass emission coefficient . When this expression is integrated the formal solution 
is 


N^(S) = N^(0)exp (- J k^p ds) +J j(S ) exp (- I k^p ds)ds 

o o s' 

The first term on the right is the decrease in radiance due to atmospheric 
extinction, and the second term is the increase due to emission and scattering. 


The source term is very difficult to evaluate for a scattering atmosphere 

since it is necessary to completely characterize the aerosol and radiating sources 

in order to perform the geometric scattering function integration . 


For a non-scattering atmosphere the emission coefficient is considerably 
simplified, and can be written simply as 


j 


v 


(S) = k (S) B (T) 

v v 


he 



When (T) is the; plank radiance function describing a black body at 
temperature T. 

In the usual atmospheric temperature case, where one does not have induced 
emission, the transmission function approach is a convenient method of 
evaluating spectral emissivity e ^ of the atmosphere. 


V s > 



i 


V s > 
V T > 


If for the moment one neglects the radiance (0) into the gas layer, and assumes 
an isothermal atmosphere 


T (S) = 1 
v 


S S 

- \ exp (- \ 

c I 


exp (- \ p ds) ds' 


-i 


= exp (-r \ p ds) 


Hence one may write 



APPENDIX B 


RELATION BETWEEN RADIATION SPECTRUM AND INTERFEROGRAM 

In this section the relationship between the spectrum of the radiation incident 
upon an interferometer system, and output of the system will be outlined. 

Incident Radiation 

The radiation incident at the instrument will have a dependence with spectral 
frequency which may be conveniently represented 

S(v, Q) = S q X T (V, Q) . 

The vector Q has been used to represent all the parameters which determine the 
detailed structure of the spectrum; things such as quantities and profiles 
of various absorbers along the optical path, temperature pressure structure along the 
path, gradient of any non-flat reflectors: (in short, anything which results 

in structure) . The results of all such parameters may be considered as a dimension- 
less transmission factor T(i),Q). The received spectral radiance S and the background 

2 

radiance S on the otherhand have dimensions of (watts/cm sterad. wavenurr.' ers] . 

o 

Defined in this manner, S is independent of the interferometer system. However, 
it is convenient to include in Q the characteristics of any spectral filter 
used in the instrument to isolate the region of interest in the spectrum. Thus 
T will include the filter transmission function. 

INTERFEROMETER OUTPUT 

The remainder of the interferometer system is characterized by its light grasp 
or etendu E (cm steradians) , its efficiency e (dimensionless) and of course the 
delay setting x (centimeters) . (The efficiency e may be used to account for 
variation in x over the interferometer aperture by including a fringe ! 
visibility term.) For a given setting of the interferometer, the radiation 
available to the detector is * 



I (x, Q) = e E 


00 

J t 1 + 


cos (2ttvx)] S(v,Q)dv 
V = 0 CO 


e E S [ 
o 


1 T(v, 
v=0 


Q)dv + | cos (2rcvx) T (v,Q)dv] 


We may define = e E S q (nominal available spectral power, 

[watts/wavenumber]), and 


W(x,Q) 



( 27 tvx) T (v ,Q)dv 


(equivalent width of optical path + interferometer system for delay setting x, 
in wavenumber) . 


Then I(x,Q) = P q (W(0,Q)] (watts). 


The equivalent widths W (x,Q) depend only on the characteristics of the 
optical path and the spectral filter, and may be evaluated independently of 
considerations of the gross parameters S q , e, and E. 

Expansion of Spectrum 


It is generally desired to operate under a range of conditions for Q. Signal 
processing will have as its object the measurement of one or more of the 
parameters comprising Q, in such a manner as to be insensitive to variation of 
the remaining parameters within their anticipated range. In developing a procedure 
for making this measurement it is useful to consider approximations for the 
variation of the spectra and interf erograms with the parameters Q. 

If nominal conditions are defined by parameters Q^, it is possible to write 
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S(v,Q) = S(v,Q ) + E(Q i - Q ) S i (v) + 


s 00 (v) + 1 <Ji S i (v) + E S ij (v) + 


Where ; if the terms of such an expression are obtained from a Taylor 
expansion, one has 


. , 83 ^, - 


3 8 . 

1 


ij ‘SQ/Qj ‘ *o 


etc. 


In the case of Q representing the concentration path products for absorber 
with absorption coefficients or, one has 

S(v,Q) = S o e" ZQ i a i 

* <= S e J ‘®io° t i [1 - (Eq.a.) + (£q.a.) 2 /2l ] 

O X JL IX 


S + I q, S, + 
oo 1 i 


This gives 


_ _ -IQ. a. 

S = S e *io i 

oo o 


and S . = S a . 

i oo x 


Expansion of Interferogram 


Because of the linear property of fourier transformation/ any expansion of the 
spectrum as a sum of components results in an interferogram which is the sum of 
the interferograms of the individual spectral components. ‘ 
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Thus, 


S(v,Q) = S oQ (v) + Iq. S i (v) + Zq ± q j (v) + 


gives an interferogram 


I(x,Q) = I OQ (x) + I q. I i (x) + E q ± qj 1^ (x) + 


00 


where I Jx) = i [1 + cos (2ttvx)] S (v)dv , 

00 K=o 00 


It is instructive to consider the case where only the linear terms are important. 

If there are m parameters q^ (i = 1, . . .m) , and if the interferograms is 

observed at n points x (K = 1, ... n) , then the observations are 

K 

I(X 1 ) = X oo (X 1 } + q i I l (x l ) + <3 W 1 m (x l ) 


Kx ) = I (x ) + q I (x ) 
n oo n 1 1 n 




m n 


If there are just m observation points, and if the I oo terms are- known, then 
the observations consititute m equations in m unknown (the q J , and may be solved 
by matrix inversion. The equations may be written 

I( V ’ r oo <V 5) - ’ , % W 

1=1 

•= E q. I .. 

, ■ ik 


which have the solutions 


m 


„rl 


’i-J , I ijt tI( V - r oo <V 2” 


- 5-1 



In the case that there are more observation points than there are variables/ the 
linear least squares technique could be applied. 




APPENDIX C , 


DETERMINATION OF WEIGHTING FUNCTIONS 


Basic Principles 

I 

In the linear approximation, it may be assumed that the interf erogram may be 
represented in terms of the interferogram obtained for some nominal set of 
conditions (nominal atmospheric parameters and nominal composition) , to which 
is added differential interf erograms according to the deviations of the conditions 
from nominal. For example if the parameters which determine the spectral 
distribution of the incident light are represented by some array of values Q, i.e. 

S = S(v, Q) , (Q = Q^/ Q 2 # Q^) 


then at least for small deviations of Q from some nominal set of values Q we may 
write 

N 

S(v, Q) = S(v, 5) + r q. (dS/dQ. ) ~ 
o i X Q 

1 ° 

= S(v,Q ) + Eq. S.(v), 
oil 

0 

I 

where q. = Q. - Q. . 

^i i x io 

Now the interferogram observed I(x,Q) is the cosine fourier transform (CFT) of S(v,Q). 

r . : 

As Fourier transformation is a linear operation, we may write to the same 
approximation 

N 

I (x,Q) = I (x,Q ) + Z q. I . (x) , 

1 

where. I^(x) are the CFT's of S^(v). 

(Non=linear effects may be included by adding terms of the form q.q.I .(x)). 

1 3 ID 

It is desired to measure one of the parameters, say. Unfortunately variation 
of the other parameters from nominal values as well as this one result in variation 
of the interferogram. The problem is to make a measurement on the interferogram 
which is sensitive only to variations in the target parameter, q^ , and not to the othe 
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Signal processing theory has dealt with the application of a linear filter, 
or weighting function H(x) to the interferogram to make the measurement. An 
additional degree of freedom may be available in the length of time D(x) (duration 
spent at each interferogram point (D(x) must be positive). The measurement made 
is then 

M(Q) = Jh(x) D (x) I (x ,Q) dx 

= M + 2 q . M . ; 
o 1 r 

where M. = H(x) D (x) I . (x) dx, 

1 f 1 

M = H (x) D (x) I (x,Q )dx. 
o o 


It is relatively simple to produce a weighting function which gives zero 

for all the M^, except for the target M ^ . It is only necessary the the range 

of x considered contains at least as many points which exhibit different 

variation (with the parameters Q) as there are parameters. This requires that the 
• ! • 
matrix B. . 

JO 


3 


J - Pi 


(X) 


1^ (x)dx 


be non singular. 


The wide, degree of freedom remaining in the choice of a filter function 
(after requiring = 0 for i ^ j) may be used to minimize the noise 
or uncertainty obtained in making the measurement. 



Representation of Noise 


The noises can be divided into four groups, according to whether they are 
random in time, or synchronous with the scanning of the interferometer delay, and 
according to whether they are basically additive or multiplicative in their 
affect on the interferometer output. The table below gives examples of sources for 

I 

the various noises. 



Random : 

Synchronous 

Additive : 

1. Detector noise, photon 
noise ; 

3. Unanticipated interf erents ; 

Multiplicative : 

2_ Random vibrations, 
scintillation; 

4. Synchronous vibration, scan 

sweep rate deviations , incorrect 
specification of weighting 
function. 


The difference in these noise sources is due to the way in which they affect 

the measurement. Any point of the inter ferogram x supplies a contribution to the 

measurement 

Am (x) = dM/dx 

*= H (x) D(x)I cj^l^(x). (M= Jam (x)dx) . 

Associated with this contribution there is an uncertainty, AN(x), which adds in an 
RMS manner with the uncertainties or noises at all other points to give the total 
uncertainty. 

2 P 2 1 

N = J(An) dx, 

For synchronous noises, AN(x) increases directly with the length of time the 
point is observed (the duration D(x>). For non-synchronous noises, AN increases 
only with the square root of D. I.E. 

J 

AN(x) «{D(x) synchronous noises 

{d^^x) non-synchronous noises. 




In the case of additive noises, AN is independent of the signal amplitude , 
while in the case of multiplicative noises N is proportional to the expected 
size of the signal at that point: 

AN « <I(x)> 

- a <qj> I^(x)) 1/2 

The notation 

r 2 (x) = E <q 2 > I 2 (x) 

will be used in representing this portion of the multiplicative noises. 

In addition to the factors outlined above, each of the contributors to A N is 
weighted by the full weighting function H(x).. The various noises will again 
add up in an RMS manner; i.e. 

2^2 
AN (x) = E AN^ (x) . 

k=l * 

The table below then shows the way in which the various noise sources contribute 
2 

to AN , (apart from constant factors which indicate the total severity of the sources 


Variation of Noise Sources AN^(x) 



Random 

Synchronous 

Additive 

1. H 2 (x)D(x) 

2 2 

3. H (x)D (x) 

Multiplica- 

2. H 2 (x)D (x) r 2 (x) 

4. H 2 (x) D 2 (x) r 2 (x) 

tive 





If the constant factors are a^ (i = as indicated in the table) , then 

we may write 


where 


AN 2 = H 2 [(a 1 + a ^ r 2 )D + (a 3 + a 4 r 2 )D 2 ] 

= H 2 D [G^x) + G 2 (x)D] ; 

G 1 (x) = [a 1 + a 4 r 2 (x) ] 

2 

G 2 (x) = [a 3 + a 4 r (x) ] . 


The total expression for the noise is then 




(x) D (x) [G 1 (x) + G 2 (x) D (x) ]dx; 


where G^ = a^ + a 2 r , 


G 2 " a 3 + a 4 r ' 
and r 2 = 2 <q 2 > i 2 ( x ). 


Optimization of M/N 


.Of the various factors which affect the signal to noise ratio M/N achieved in making 
measurement, only the weighting function H(x) and the duration function D (x) may 
be varied in an attempt to optimize the ratio once the basic system has been 
specified. 



It is fairly straightforward to find a weighting function H(x) which is optimum 
for a specified duration function D(x). This will be outlined now. 

It is desired to find a weighting function H(x) which gives 


~ J H (x) D (x) I^(x)dx 

0 i*l 

" 1 i=l 


This will give a measurement 


M = J H D q. I. 


dx 


H is to be chosen so as to minimize N 


N 2 « J\h (x)D (x) ] 2 G (x)dx 


where G(x) = G^/D + G 2 


Defining H' (x) = H(x)D(x), the problem is to minimiz;. 

Jh'(x) 2 G(x)dx = N 2 (1) 

S.t. Jh' ( xji^xjdx = 6 il . (2) 


The equations (2) may be multiplied by Lagrange multipliers nu and added to equation 
(1) , to require 

Jh ' 2 Gdx + Eiru J*H' I^dx = minimum . 


* 
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arying any particular H (x) will result in zero variation of this expression when 
H' is optimized: i.e. . 

3H' (x ) ^ ^ + ^ ' 1^ dx] = 0 for all x^, 

or 2H'(x o ) G(x q ) + 2m i I i (x o ) =0. 


This gives 


H' (x ) = Em! I (x )/G(x ) 
o 1 1 o o 


where m! = -m./2. 

i i 


The nu may be determined by substituting for H’ in the equations (2) and solving: 


I' 


H* (x) I. (x) dx = 6,, 
3 il 


“ £m f(I I /G) dx 

- -iV'V 

where A_ = Jdd./G) dx 


This sfet of linear eqv .tions has the solution 


' -1 

m . = I A . . 6 . , . 

1 j ]1 


where A * is the inverse of matrix A. 


Thus we may write 


H' (x) = HD = E A - * I i (x)/G(x) 


or H(x) = E A"^ I j .(x)/[G 1 + DG ]x. 


* 
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This II (x) can be shown to be the same one which would be used to find those 

values of q which give the best fit to an observed interf erogram I (x) , i.e. 

those q which minimize 

|[I(x) - Eq. I . (x) ] 2 /N 2 (x) dx, 

V 1 1 

when N 2 (x) = G^xJ/D + G 2 (x) = G. 

The following points should be noted. 

- If a given function H(x) is used, function H^(x) = CH(x) may be used to give a 
larger measurement = cM; but the signal to noise ratio remains unchanged 
(M 1 /N 1 = Mj/N^ . 

- If there are synchronous noises only, (G^ << G 2 > , then any combination of H and 
D which give the same product HD will give the same signal to noise ratio. In- 
creasing observation time in such a case does not improve the signal to noise 


ratio. 
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APPENDIX D 
COMPUTER MODELLING 


Digitize \ 
Migeotte \ 
Spectrum / 


Program REWRT 
Digit to Create 
Wave No 
Transmi ttance 
Flow 


HC1 Tape 
File From 
York. U. 


Program RDM AG 
to Create Wave No 
Trans. File 


Disk Files 
[Migeotte & HClj 



1 

Program COSPEC 
To Model Corre- 
lation Spectrometer 
Signal 





Printed Listing 
COSPEC Signal 



Program I NT ERF 
To Compute Fourier 
Transforms of -IP 
Files 


Disk Files 
.Fourier Transforms 


Program RUNEM 
Create Weighting 
For & Model Signal 
Of Interferometer 


Printed Listing 
Correlation 
Interferometer 
Relative Signal 











DATA PILES 


Cl) Atmospheric spectra profiles were obtained from The Solar Spectrum 

From 2.8 to 23. 7y Part 1, Photometric Atlas (Migeotte, Nevin, Swenson) 


These profiles were digitized onto aards and placed onto disk files'. 

The disk files were run through programme RF.WRT, which reversed the order of 
the data to allow conversion from wavelength to wavenumber; and programme 
DIGIT which eliminated digitizing errors by linearizing and rotating 
X-axis, splined Y-axis data to equal X axis intervals (.025 cm ) 
and output a disk fit of wavenumber versus transmittance. 

Programme NORML normalized data to maximum Y = 1.0 and created 
a binary disk file across the wavenumber range of interest (2700 - 2910 cm” 1 ) . 

(2) A tape file of wavelength versus transmittance for HC1 gas of various 

concentrations was obtained from a radiative transport computation done 
at York University. This tape file was run through programme RDMAG to 
convert the data from IBM to PDP equipment compatible mode, convert 
from wavelength to wavenumber, and spline to equal frequency intervals 
(.025 cm ) . The new files were output onto disk over frequency range 
2700 to 2910 cm" 1 . 


This procedure was repeated for each of three of the original tape files. 
Instrument Modelling 

(1) Non-dispersive Analyzer (NDA ) . Programme ANALY-+ calculates the 

signal through the NDA due to HC1 gas in the atmosphere by calculating the 
zero-gas correlating factor C = E ^M * IF* HCl cg ^ ^ / E(M*IF) 

where M is Migeotte solar spectrum transmittance. 


IF is interference filter transmittance, and this was chosen to be the function 
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2 3 —1 — 1 
(l-(F/Fo) ) where F is frequency 1 cm F _< 200 cm and Fo = 100 cm 

This function approximates what may be expected for about a 90 cm half • 

width filter, and has a finite extent for ease in computing. 

HCl cell is the HC1 gas cell transmittance of a cell containing 0.1 
atmospheres of HC1, and 1 cm length at 295°K. 

The signal was computed as follows 

£ CM * IF * HC1 ^ * HCl )~C*£(M*IF* HCl 1 , ) 

atm cell atm 

Where HC l atm i s HC1 9 as the atmosphere with concentrations of 

1 ppm-m or 100 ppm-m. The programme also calculates the signal 

using the square root of the transmittance of Migeottes spectra as a measure 

of changes in signal due to different atmospheric amounts of interferent 

species. 

(2) Correlation Spectrometer 

The programme SPECT calculates signal resulting from correlation 
of two exit slit masks. 

The programme sets up two sets of elements in the exit slit, MSKl and 
MSK 2, at specified wavenumber locations which are on and off the HCl gas 

i ' 

peak respectively. 

The signal is then calculated for the gas and no-gas conditions as 

Sgas =[E(M * HCl * IF). MSKl] - 21 CM I HCl * IF). MSK 2 ) 

which is resultant signal due to Migeotte solar spectrum (M) , HCL 
atmospheric gas, and interference filter (IF) as before and 

S 

no gas 


£(M * IF). MSK 1 


] -[ 


£ (M * IF) . MSK 2 


Is signal with no HCl gas present. 



The output signal of COSPEC due to HC1 gas alone is then 

CS - S ) 

gas no gas 

This output signal is calculated with 1 ppm-m and 100 ppm-m of 

HC1 gas and also for a solar spectrum J M Ci.e. change in atmosphere) . 

The mask elements are the triangular spectrometer response functions 
resulting from equal width entrance and exit slits (1.73 cm ) 

(3) Interferometer 

The programme calculates the Fourier Transform of . the input transmittance files 
which yields both the amplitude and the phase of the interferogram 
output signal. 

The programme applies the same interference filter over the input 

spectral range 2705 - 2905 cm -1 as used for the previous sensor calculations, 

(1 - (F/Fo) 2 ) 3 . 

This filtered data is then interpolated using the smoothing function 

( sin ( 2tt f c n ) j ( sin 2tt F^ n 
2?r f c n / \ 2tt Fn n 



Where f is cut-off frequency 

P i s the apodizing frequency 
n 

n is an index of frequency location 1 _< n _< N 

This reduces the input data 1024 points. The fast Fourier Transform 
(FFT) is then applied to the 1024 points using the Cooley-Tuley 
algorithm, and files written into disk. 


The FFT files are used by programme RUNEM to create a weighting function 



which, was used as the correlation function for the test interferograms . 


The following library set of interferograms were computed: 


51 = E(Migeotte * Interference filter) 

52 = e/m * IF * HC1 (1 ppm-m)) 

53 = zlrvT * IF * HC1 (1 ppm-m)) 

54 = e( tir * if) 

55 = 2^M * IF * HC1 (100 ppm-m)) 

56 = E (nr * IF * HC1 (100 ppm-m)) 

A few samples of computed interferogram amplitudes are shown in Figure Dl. 

The weighting function is created using the library interferograms SI, S3 and S4. 

Table Dl shows the results of the application of the derived weighting functions to 
a range of interferogram segments. The uppermost number in each case is the 
deduced HC1 burden deduced from S3 relative to S2 which had 1 ppm meter of HC1. 

The best interferogram segment should record 1.00 ppm-m. -As seen from the 
table several regions allow a measurement of the expected burden (within 
a few percentage points) . Furthermore the second and third numbers in each 
set is the computed burden for 100 ppm-m HC1 burden. The number returned 
for a perfect measurement with a linear curve of growth would be 100 in every 
case. As can be seen the curve of growth is not quite linear. Further 
the third number is deduced for a different interferent spectrum from 
the second number. The table dramatically illustrates the interferent 
rejection capability of correlation interferometry analysis. 
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RELATIVE HC1 MEASUREMENT OUTPUT 


HCl RELATIVE S/N 



RUNEM MEASUREMENTS ON FFT ' S S3, S5, S6, RESPECTIVELY FOR VARYING DELAY INTERVALS OF CORRELATION 

(RELATIVE TO S2 = 1.0) 


TABLE Dl- 


Continued 



DELAY (cm) 


FIGURE Dl. COMPUTED 


INTERFEROGI 


__ Migeotte x IF " (*S I ) 

Migeotte x IF x IOO ppmm HCL (S5) 

__ 100 ppmm HCL x IF 



SEGMENTS 



